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ABSTRACT 


The  analytical  and  experimental  results  for  traveling 
bending  waves  in  a  uniform  beam  with  tuned  periodic  canti¬ 
lever  ribs  are  presented  in  this  report.  Two  analytical 
models  were  developed:  one,  is  for  an  infinite  array  of 
periodic  sections  and  the  other  is  for  a  finite  array  of 
periodic  sections  in  an  infinite  beam.  Several  types  of 
periodic  sections  were  considered.  The  bending  moment 
spectrum  was  measured  for  a  traveling  wave  in  a  uniform 
beam  caused  by  impact.  The  insertion  loss  as  a  function 
of  frequency  was  measured  for  such  a  beam  with  1,3, 5, 7 
and  15  sets  of  symmetric  sections.  There  is  good  agree¬ 
ment  between  the  computed  bending  moment  spectra  and  the 
measured  ones  for  all  sets  of  cantilevers  considered.  The 
width  of  attenuation  bands  increased  with  an  increasing 
number  of  cantilever  sets,  and  the  spectrum  for  15  canti¬ 
lever  sets  was  a  close  approximation  to  that  of  the 
computed  spectrum  for  an  infinite  number. 
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1.  INTRODUCTION 

The  theory  of  wave  propagation  in  periodic  structures  has 
been  developed  by  solid  state  physicists  and  electrical 
engineers  for  some  time  [1,2}.  The  application  of  this 
theory  to  predicting  vibration  responses  of  periodic 
engineering  structures  due  to  broad  band  pressure 

fluctuation  is  relatively  recent  and  has  been  applied, 
primarily,  to  the  analysis  of  periodic  structures  as  band 
pass  filters  for  bending  waves  £3, 4, 5, 6].  The  analysis  for 
an  infinite  ribbed  bar  has  been  published  with  some  limited 
confirming  evidence  C7l-  The  objective  of  this  study  is  to 
get  extensive  quantitative  results,  both  analytical  and 
experimental,  for  a  beam  with  tuned  periodic  cantilever 
ribs.  Two  analytical  models  are  presented  for  such  a 
periodic  structure  in  which  the  periodic  sections  consist 
of  constant  cross-section  rectangular  beam  elements.  One  of 
the  models  is  for  an  infinite  array  of  periodic  sections 
and  the  other  is  for  a  finite  array  of  periodic  sections  in 
an  infinite  beam.  The  attenuation  in  the  infinite  model  is 
either  zero  or  infinite?  the  stop  and  pass  bands  are 
conveniently  found  by  an  eigenvalue  technique.  The  actual 
value  of  the  attenuation  must  be  found  for  a  finite  model 
and  so  an  impedance  technique  is  used.  Three  types  of 
periodic  sections  are  considered,  symmetric  (two  sided 
cantilevers) ,  asymmetric  (  one  sided  )  and  antisymmetric 
(  alternately  on  each  side  ) .  The  impedance  model  is 


sufficiently  versatile  to  predict  the  response  of  ordered 
as  well  as  disordered  periodic  sections  with  no  limitations 
on  the  number  of  such  sections.  This  report  gives 
experimental  and  analytical  results  for  an  ordered  periodic 
beam  with  symmetric  sections  of  1,  3,  5,  7  and  15  sets  of 
cantilevers. 

In  an  earlier  report  the  analytical  model  for  the 

infinite  periodic  beam  with  symmetric  sections  was 
developed  and  compared  with  measured  results.  This  work  is 
briefly  reviewed  in  Section  2.3  of  this  report. 

In  this  report.  Section  2  gives  the  analysis.  Section  3 
describes  the  experimental  method.  Section  4  contains  the 
results  and  Section  5  discusses  the  comparison  of  the 
experimental  and  analytical  results. 

2.  ANALYSIS 
2 . 1  GEOMETRY 

Our  periodic  structure  consists  of  an  array  of  priodic 
sections.  Elements  of  the  periodic  section  are  uniform 
beams  of  rectangular  cross-section.  Three  section  types  are 
considered; 


a)  Symmetric:  Sections  of  this  type  consist  of  a  uniform 
beam  element  with  a  double  cantilever  element  attached  to 
it  at  its  right  end,  figure  A.l. 


b)  Asymmetric:  Sections  of  this  type  consist  of  a  uniform 
beam  element  with  a  cantilever  attached  to  one  side  at  the 


right  end,  figure  A. 2. 

c)  Antisymmetric:  Sections  of  this  type  consist  of  a 
uniform  beam  element  with  cantilever  elements  attached 
alternately  on  each  side,  one  at  the  center  and  the  other 
at  the  right  end,  figure  A. 3. 

2.2  EQUATIONS  OF  MOTION 

In  first  order  theory,  the  differential  equation  for  the 
propagation  of  an  extensional  wave  in  a  uniform  beam  is: 

zfu  1xu. 

ea  T71  +CA  at1-  =  °  (2.2.D 


and  that  for  a  bending  wave  (  Bernoulli-Euler  )  is: 


VS  VS  _ 
D  at1-  “ 


(2.2.2) 


If  we  assume  a  sinusoidal  wave  of  frequency  u>  travelling  to 
the  right 

u(x,t)  =  u  ( x )  exp  (-ju>t) 

(2.2.3) 

y  (x,  t)  =  v  ( x )  exp  (- j  lOt) 

where  j  =  V  ~1 
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the  general  solutions  of  equations  (2.2.1)  and  (2.2.2)  for 
a  uniform  beam  are 


u(x) 

v  (x) 


EjCOS  Jfx  +  E  sin  jfx 
A  jCosh  p  x  +  A^sinh^x  +  A^cos 


(2.2.4) 


x  +  A  sin 

4 


f* 


where 


V  =  <*>/  cr 
cr®  V  E/f  is 

j9l+=  12/>uf/ 


the  rod  phase  velocity 
Et 


t  *  beam  thickness  for  a  rectangular  cross-section 
D  =  El  is  the  bending  stiffness*. 


Aj ,  A^,  Aj,  A+,  E4  and  E^  are  constants  to  be  determined 
from  the  end  boundary  conditions. 


We  can  analyse  a  travelling  wave  in  a  periodic  structure  by 
using  the  nodal  boundary  conditions  at  the  ends  of  the 
periodic  sections. 


2.3  TRANSFER  MATRIX 

We  define  a  vector  composed  of  six  components  of 


*  the  anticlastic  curvature  restraint  is  not  effective  for  a 
thin  rectangular  cross-section  at  small  bending  strains  and 
so  the  1-t>  term  is  ommited. 


displacements  and  forces  at  a  node  i 


Iv 


zr 


u , 


EA 


0c /p 

(_-F;  /El  f'a  J 

where  the  components  are  defined  in  appendix  E  and  figure 

2.2.1. 


(2.3.1) 


The  elements  of  the  vector  Z^  are  multiplied  by  suitable 
factors  to  give  each  the  dimension  of  length.  We  can  find 
the  above  vector  at  the  end  of  a  section  in  terms  of  the 
that  at  the  other  (figure  A.l,  A. 2  or  A. 3): 

zi  "  fH Z^.,  (2.3.2) 

[H,J-  is  a  (6x6)  matrix  and  is  the  transfer  matrix  for  the 
section  i,  between  nodes  (i-1)  and  i  of  section  type  m 
(m=l,2,3).  We  partition  the  matrix  to  get 


where  H,m  (2x2)  and  H4m  (4x4)  are  due  to  the  axial  and 
bending  modes  alone  and  H2W  (2x4)  and  H3w  (4x2)  are  the 
matrices  coupling  the  two  modes  of  vibration.  The  transfer 
matrices  for  each  type  of  section  are  given  in 

appendix  A. 

2.4  EIGENVALUE  MODEL 

The  periodic  beam  is  considered  to  have  an  infinite  array 
of  identical  periodic  sections.  If  we  have  a  wave  of 
frequency  co  traveling  in  an  undistorted  fashion,  it  must 
attenuate  by  the  same  fraction  in  each  section.  Thus  we  may 
write 

z;  =  “Xyn  ZC_,  (2.4.1) 

where  A*,  is  the  ratio  between  successive  periodic  repeats. 
From  equation  (2.3.2)  we  have 

Zi  =  t«m]  Z  i-i  =  **  ZH  (2.4.2) 

This  is  an  eigenvalue  problem  where  [Hm’]  is  a  function  of 
frequency*.  Since  all  of  the  sections  are  identical,  the 

*  submatrices  G  and  F  of  [Hm].  were  incorrectly  reported  in 
reference  [8] .  The  corrected  version  is  given  in  equation 
(E.2)  and  (E.7),  appendix  E. 
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subscript  i  of  the  transfer  matrix  is  dropped.  There  are 
six  roots  to  the  eigenvalue  problem  which  are,  in  general, 
complex.  It  can  be  shown  that  the  absolute  values  are  equal 
to  or  less  than  one.  If  the  value  is  one,  the  magnitude  of 
the  travelling  wave  is  unattenuated;  if  it  is  less  than 
one,  it  must  eventually  die  out.  The  stop  and  pass  bands 
calculated  for  h/L  =  0.48  are  shown  in  figure  4.1.1  and  the 
predicted  results  are  compared  with  measured  ones  in  figure 
4.1.2  which  is  discussed  in  section  4. 

2.5  IMPEDANCE  MODEL 

The  structure  is  modeled  in  three  parts.  The  middle  part 
of  the  structure  consists  of  a  finite  array  of  a  number  of 
periodic  sections;  the  two  end  parts  are  semi-infinite 
uniform  beams  attached  to  the  middle  part.  Since  the 
experimental  portion  of  this  project  considers  the  response 
of  the  beam  to  the  impulsive  force  applied  transversely 
near  the  center,  the  Sommerfeld  condition  requires  us  to 
consider  only  waves  travelling  to  the  right  at  the  right 
end  and  to  the  left  at  the  left  end. 

2.6  IMPEDANCE  MATRIX  and  INSERTION  LOSS 

The  impedance  matrix  Q*  relates  the  force  vector,  N^,  at  a 
node  to  the  displacement  vector,  W^,  at  that  node  (figure 
2.6.1).  The  sign  used  as  the  superscript  on  the  impedance 
matrix  signifies  the  direction  in  which  one  looks. 
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Where 


N; 


Mj  /El  |S4 
-  -Ft  /El  £3  > 

PL  /  IS  EA 


W; 


We  now  rearrange  the  vector  ,  for  convenience,  as 


(2.6.1a) 


The  impedance  matrix  at  node  n  which  is  the  first  node  of 
the  right  hand  semi-infinite  uniform  beam  (figure  2.6.1)  is 
(appendix  B  equation  (B.l)): 


-3 

- (1- j) 
0 


Nn  =  Q 


*'V  r  . 


-d+j) 

j 

0 


r,  «r, 
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;  j=  4~ i 


(2.6.1b) 


A  ^  ' 


The  transfer  impedance  matrix,  [R*r>]w  gives  the  force  vector 
at  the  n’th  node  in  terms  of  the  displacement  vector  at  the 
0 ' th  node : 


Nn  =  Jm  We  (2.6.2) 

It  is  shown  in  appendix  C  that  [R*n]  can  be  computed  from 
,  the  inverse  of  the  transfer  matrix: 

-I  H 

[  {*"»']*,  +  Qn  1  [  1  •  •  •  • 

•••[W2  +  iK™]xQV}  [1M,  +{Kwl}t  Q  +  ]'1  (2.6.3) 

A  bending  wave  traveling  to  the  right  in  a  bar  will  split 
into  a  reflected  wave  and  a  transmitted  wave  at  a  change  of 
cross  section  which  occurs  at  a  node.  Therefore  at  the 
zero' th  node 


( 


or 


W1  "  W*  +  WT  "  W*  +  W0 


(2.6.4) 


In  the  same  way  for  the  forces 


or 


NI  -  NR  +  *T  "  N*  +  NC 


(2.6.5) 


* 

y  * . 


The  relations  between  force  and  displacement  vectors  al 
node  0  are 

toll"  .  V  [oj] 


N  =  Q+  W  =  Q+  W. 
t  or  o  o 


(2.6.6) 


where 


Mfc/EI  pZ 

'  Nk  *  ^-pn/EI^S 
P*  /JTEA 


(k»I,R,T) 


[.Q^]  is  obtained  from  equation  (D.5),  appendix  D. 


Q*  is  obtained  from  equation  (C.4) ,  Appendix  C  for  i=0. 
Since  the  reflected  wave  at  node  zero  propagates  through 
the  uniform  semi-infinite  beam  in  negative  x  direction, 
is  obtained  in  the  same  way  as  Q*is  obtained  in  appendix  B 
but  with  opposite  sign  of  x. 

Using  equations  (2.6.6)  and  (2.6.4),  equation  (2.6.5)  can 
be  written 


«,  -  -ton,  [o;f + 

m  {r  -  (i  -  [q'.jjq:]"}  no 

ni  -  { i  -[otijot]"  K «. 

Ni  *  s  w. 

where  S  -  {l  - [0^r [Q.* j‘' j '{l  -  [o;jH j 0* 

So  that  WA  =  s"1  N._  (2.6.7) 

o  I 

S  is  a  (3x3)  matrix,  given  in  appendix  D.  Substituting 


equation  (2.6.7)  into  equation  (2.6.2)  we  obtain, 


N  =  J.N-  (2.6.8) 

M  TTl  X 

where  J_  =  [R*  }  S  1 
m  1  o*Jyr) 

Since  bending  strains  were  measured  in  our  experiment  we 
consider  the  bending  moment  element  of  Nh  in  equation 
(2.6.8) : 

Mt./Eip2  =  Jm(l,l)  M,/EI^+  Jw  (1,2)  .-F3 /EI^3  +  Jw(l,3)  P^/tfEA 

Substituting  the  values  of  Mx  ,  and  in  terms  of  the 
amplitudes  of  the  incident  axial  and  bending  waves 
(equation  (D.5)  we  find  that 

Mn/EIf3*=  [-Jm(l,l)  +  jJm(l,2)]  V  +  j J-n, (1,3)  U,  (2.6.9) 

For  the  beam  with  symmetric  sections,  J  (1,3)“0  so  that  the 
bending  mode  is  not  coupled  with  the  axial  mode. 


The  dB  Insertion  Loss  (IL)  due  to  the  presence  of  the 
periodic  section  is  = 


IL(w)  =  20  logj— |=-20  log(e  ||  Jm(l,l)  -  jJ„,(l,2)||  (2.6.10) 


We  have  assumed  that  the  damping  dissipation  in  the  uniform 
beam  is  negligible  so  that  the  spectrum  magnitude  is  the 


same  everywhere.  Spectrum  measurements  confirm  this 
assumption.  Thus  the  transmitted  spectrum,  Sp ,  can  be  found 
from  the  incident  spectrum, 

Sp(w)  =  Su(i»)  -  IL(*)  (2.6.11) 

where  Su (w)  is  the  uniform  beam  spectrum. 

3.1  INST  RUMENT AT I ON 

With  one  exception,  the  instrumentation  is  the  same  as  that 
described  in  reference  [8],  figure  3.1.1.  Instead  of  using 
the  HIPlot  X-Y  recorder  to  make  hard  copies  of  the  signals 
stored  in  the  spectrum  analyzer  a  HP-85  computer  and  a  HP 
7  225B  plotter  were  used.  This  system  is  more  flexib.’e  for 
getting  a  hard  copy  of  any  stored  signal. 

3.2  BEAM  and  CANTILEVERS 

A  uniform  steel  beam,  0.250  inch  thick  by  2.00  inch  wide 
and  36  feet  long  was  hung  in  a  horizontal  position  by  0.030 
inch  piano  wires  fastened  at  4  feet  intervals  to  the  thin 
edge  by  welding  to  3/4  inch,  6-32  screws  in  drilled  and 
tapped  holes;  the  screws  were  locked  by  locknuts.  The  wires 
were  fastened  to  an  overhead  pipe  4  feet  above  and  adjusted 
by  turnbuckles  to  equal  tension  by  tuning  to  the  same  note. 
The  impact  point  was  16  feet  from  one  end.  The  travelling 
waves  were  measured  at  point  10,  10  inches  from  the  impact 
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point  and  point  82,  82  inches  from  the  impact  point,  both 
on  the  20  feet  end  (figure  3.2.1). 

After  the  uniform  beam  measurements,  15  pairs  of 
cantilevers  each  2.35  inch  long  were  fastened  to  the  beam 
at  5  inch  intervals  starting  one  inch  from  point  10;  the 
measured  resonances  of  the  cantilevers  mounted  in  a  short 
length  of  beam  lay  between  1332  and  1386  Hz.  The  resonant 
frequency  of  a  2.35  inch  cantilever  calculated  by  using  the 
measured  phase  velocity  is  1394  Hz;  if  we  assume  the  length 
(arbitrarily)  to  be  2.40  inch  to  account  for  the  elasticity 
of  the  base,  we  get  1340  Hz.  The  frequencies  of  waves  with 
wavelengths  10  and  5  inches  (5  inches  is  the  distance 
between  two  consecutive  cantilevers)  are  respectively  863 
and  3452  Hz. 

3.3  EXPERIMENTAL  TECHNIQUE 
3.3.1  FREQUENCY  SPAN 

Since  the  bending  wave  is  dispersive  with  its  phase 
velocity  proportional  to  the  square  root  of  its  frequency, 
the  low  frequency  components  of  the  wave  take  longer  to 
travel  from  the  origin  of  the  wave  at  the  point  of  impact 
to  a  measurement  point  than  do  those  of  higher  frequencies. 
With  the  finite  beam  that  was  used  in  our  experiment,  the 
spectrum  of  an  infinite  beam  can  be  obtained  if  the 
spectrum  is  recorded  before  any  reflection  from  the  free 


end  of  the  finite  beam  alters  the  spectrum;  that  was 
achieved  with  signal  travel  times  appropriate  to  the 
frequency  span  settings  of  the  spectrum  analyzer.  For  a 
frequency  span  of  10  kHz  (12.5  msec  time  signal)  with  6.5 
msec  pretrigger  time,  the  total  travel  time  (12.5-6.5  or 
6.0  msec)  to  reach  the  measurement  point  82  (figure  3.2.1) 
was  not  enough  for  waves  of  frequency  10.4  kHz  and  lower  to 
return  after  reflecting  from  the  free  end  A  of  the  beam. 
Also  it  was  not  enough  for  frequency  500  Hz  and  lower  to 
reach  the  point  82  directly  from  the  impact  point.  So  at 
this  frequency  span  setting  the  spectrum  between  500  Hz  and 
10kHz  was  essentially  the  same  as  that  of  an  infinite  beam. 
At  the  5  kHz  and  2.5  kHz  frequency  span  settings  with  25 
and  50  msec  time  signals  and  9  and  17.5  msec  delay  times, 
travel  time  for  the  wave  to  reach  point  82  was  respectively 
(25  -  9  or  16msec)  and  (50  -  17.5  or  32.5  msec) .  At  these 
two  settings,  using  the  same  reasoning,  the  spectra  were 
good  from  63  Hz  to  1300  Hz  and  from  15  Hz  to  350  Hz 
respectively.  Thus,  to  get  the  spectrum  for  the  infinite 
beam  for  the  whole  frequency  range  of  10  kHz,  spectrum  data 
for  frequencies  15  -  320  Hz  and  320  -  1480  Hz  were  taken 
from  the  2.5  kHz  and  5  kHz  settings  respectively  and  the 
rest  of  the  spectrum  data  were  taken  from  the  10  kHz  span 
setting.  These  three  data  sets  were  then  combined  with 
appropriate  band  width  corrections. 


3.3.2 


UNIFORM  BEAM 


The  uniform  beam  spectrum  was  measured  at  point  82.  All  the 
measurements  were  taken  with  eight  averages  of  the  r.m.s. 
spectrum  with  gain  settings  at  10  mV  max.,  and  A.C. 
coupling.  The  bending  moment  spectrum  is  shown  in  figure 
3.3.1. 

3.3.3  BEAM  WITH  CANTILEVERS 

The  settings  of  the  spectrum  analyzer  and  the  pretrigger 
time  were  kept  the  same  as  in  the  uniform  beam  case.  The 
amplitude  of  the  impact  was  also  kept  the  same  to  make  sure 
that  the  incident  wave  was  the  same  as  for  the  uniform 
beam.  Since  the  spectrum  measured  at  any  point  beyond  the 
cantilevers  along  the  uniform  (undamped)  beam  remains  the 
same,  all  of  the  measurements  of  the  spectrum  were  taken  at 
point  82. 


4 . 1  IMPEDANCE  MODEL 

Equation  (2.6.10)  was  used  to  calculate  the  dB  Insertion 
Loss  (IL)  as  a  function  of  frequency  for  the  symmetric 
periodic  case  for  1,  3,  5,  7  and  15  pairs  of  cantilevers. 
Since  the  incident  wave  was  assumed  to  be  same  with  and 
without  the  cantilevers  the  computed  spectrum  for  the 
periodic  beam  was  calculated  from  equation  (2.6.11)  for 
each  number  of  sets.  The  computed  spectra  are  compared  with 
the  experimental  ones  in  figures  4.2.2a,  4.2.2b,  4.2.2c, 


4. 2. 2d,  and  4.2.2e.  In  some  stop  bands  the  computed 
spectrum  values  were  less  than  -120  dBV  and  were  below  the 
bottom  of  the  figures. 

4.2  EIGENVALUE  MODEL 

The  pass  bands  and  stop  bands  for  an  infinite  periodic  beam 
calculated  by  the  eigenvalue  method  reported  in  £8]  are 
shown  in  figure  4.1.1.  In  this  model  of  an  infinite  number 
of  periodic  sections,  no  attenuation  occurs  in  the  pass 
bands  and  no  transmission  occurs  in  the  stop  bands.  The 
spectrum  computed  in  this  way  is  compared  with  the  measured 
spectrum  of  the  periodic  beam  with  fifteen  sets  of 
cantilevers  in  figure  4.1.2. 

5.  DISCUSSIONS 

The  agreement  between  the  spectrum  computed  by  the 
impedance  method  and  the  measured  spectrum  is  gratifying 
for  all  five  cantilever  sets  considered.  The  measured 
attenuation  between  6600  Hz  and  7300  Hz  for  15  sets  of 
cantilevers  was  much  greater  than  the  calculated  values. 
Some  indication  of  this  can  also  be  seen  in  the  other 
figures.  The  reason  for  this  discrepancy  is  not  known; 
material  damping  is  much  too  small  to  account  for  it. 

Attenuation  in  the  stop  bands  and  the  width  of  the  stop 
bands  increased  with  increasing  numbers  of  cantilever  sets. 
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At  fifteen  cantilever  sets  they  approximated  the  computed 
results  for  an  infinite  number.  The  values  of  the  computed 
spectrum  in  most  of  the  stop  bands,  especially  for 
cantilever  pairs  five  or  more,  was  -250  dBV  or  below, 
whereas  the  experimental  attenuation  was  always  less  than 
35  dBV.  This  is  probably  due  in  part  to  the  12  bit  A/D 
converter  (72  dB  full  scale)  of  the  spectrum  analyzer  and 
in  part  to  leakage  caused  by  the  DFT . 


APPENDIX  A 


Transfer  Matrices 

Considering  the  axial  wave  solution  of  equation  (2.2.4)  for 
a  uniform  beam  piece  of  length  L,  we  obtain  (appendix  E) 

z,-t  —  ®  ^  (A.l) 

where  B  is  given  by  equation  (E.l)  and  Z(  is  defined  in 
equation  (2.3.1).  From  the  bending  wave  solution,  we  get 
(appendix  E) 


Z2i  *  G  (A. 2) 

where  G  is  given  by  equation  (E.2). 

The  process  of  getting  the  transfer  matrices  for  any  type 
of  periodic  sections  is  similar  to  that  given  in  section 
2.2  of  [8]  with  appropriate  equations  of  equilibrium  at  a 
node. 

a)  Symmetric  section  (m=l) 


The  equations  of  equilibrium  at  node  i  with  the  sign 
conventions  shown  in  figure  (A.l)  are 


p.t  - 

Pzi  + 

P»L 

P.t  + 

Fzi  - 

Fai 

Following  the  same  procedure  as  described  in  section  2.2  of 
f8]  we  obtain  the  transfer  matrix  £h.J.  for  this  type  of 

i 

section  as  (equation  E.5): 


Mi  * 


h„  ; 

i 

H2, 

B  +  D, 

1  0 

i 

i 

—  t-  — - 

j- 

i 

H„  1 

h4. 

o 

i 

i 

!  G  +  F 

• 

t 

The  matrices  Dt and  F  are  given  in  equations  (E.6)  and 
(E.7).  Here  because  of  symmetry  the  coupling  matrices 
HZ)  and  Hj,  are  null  matrices,  since  the  coupling  effect 
due  to  one  cantilever  at  one  side  is  nullified  by  the 
cantilever  on  the  other  side  (figure  A.l). 

b)  Asymmetric  section  (m=2) 

The  equations  of  equilibrium  at  node  i  are  (figure  A. 2) 


and  the  transfer  matrix  [Hz^  is 


H1Z  |  »22 

B  +  D2  j  c2 

telt  - 

1  ___ 

| 

=  _ 1 _ 

1 

HJZ!  hU2 

E2  i  G  +  1/2  F 

matrices  C2 ,  D^and  E  z  are  given  in  equations  (E.8),  (E.9) 

and  (E. 10) . 

c)  Antisymmetric  section  (m=3) 

A  section  of  this  type  can  be  considered  to  be  composed  of 

two  sections  of  asymmetric  type  where  one  is  rotated  180° 

with  respect  to  the  other  about  an  axis  along  the  length  of 

the  section  (figure  A. 3).  The  transfer  matrix  for  part  I 
l 

([Hj.)  is  the  same  as  that  of  the  asymmetric  section: 

[H,]--  [H,]; 

and  that  for  part  II  {  fH3J^ )  r  using  the  equations  of 
equilibrium 


APPENDIX  B 


Impedance  Matrix  at  the  n’th  node 

In  the  semi-infinite  uniform  beam  extension  of  the  main 
beam,  the  solution  for  the  amplitudes  of  the  time  dependent 
part  of  the  bending  and  axial  modes  of  vibration  traveling 
to  the  right  are  (figure  B.l): 

v  (x)  =  B  j  exp  {-  £  x)  +  B2  exp  (- j  p  x) 
u  (x)  =  B3  exp  (- j  K  x) 

The  above  solutions  are  completely  determined  by  the  three 
displacement  boundary  conditions  (vn  ,  Q„/p  and  un  )  at  the 
n'th  node.  We  obtain  the  components  of  the  force  vector 
(Mn/EI0z  ,  -F„ /Eip3  and  P„ /K  EA)  at  the  n’th  node  by 

differentiating  the  solutions  and  using  the  constitutive 
relations 

Mn/EIp2=  v  ^'  ( x )  I  ,  -Fn/Eip3=  vw(x)l 

X-0  'x*0 

and  Pn/y  EA  =  u'(x)| 

*X  r  0 

Since  Nn  =  Q*  Wn 

the  impedance  matrix,  Q*,  at  the  n'th  node  is 


APPENDIX  C 


C.l  Transfer  Impedance  Matrix  [Ro„] 

■'  7T1 

From  equation  (2.3.2): 


Wh 

"0 

t  Ni-i> 

•  -  r« - 

u 

where  [Km];  -  [h*,]. 

[h„J;  is  the  matrix  £hw]^  rearranged 


(C.l) 


A 

in  the  same  way  as  Zj  . 


Since  is  the  transfer  matrix  from  the  i'th  node  to  the 
(i-l)'th  node,  it  is  most  easily  found  by  changing  the  sign 
of  x  which  makes  and  «<2  negative  in  the  equations  (E.l), 
(E.2),  (E.6) ,  (E.l),  (E.9)  and  (E.10)  of  appendix  E. 


After  partitioning  the  matrix  [K^-].  ,  into  four  parts 


W 


i-i 


we  have 


Equations  (C.2)  are  valid  for  any  section  1.  Starting  wit 
the  n'th  section,  we  have: 


«V,  -  lK*»LWn  + 


Substituting  for  Nn  from  equation  (B.l)  and  inverting  the 
first  equation  to  find  Wn 


w„  -  Q t,  ]  ", 


(C.  3) 


then 


»-• 


where  Q*,  =  +  iK»«!n  O  [!*„,)„  +  Q+„] 


is  the  impedance  matrix  at  (n-l)’th  node  and  can  be 
determined  since  all  the  terms  on  the  right  hand  side  are 
known. 


In  a  similar  fashion  for  any  node  i  we  find  that 


<  -  K„};  ol,]  [(K +  <c-«> 


and  w;  =  [{"«,};  +  iK-">h  °t]" 


-26- 


in  terms 


Applying  equation  (C.4)  successively,  we  find  Wn 

of  W  .  Then  we  find  the  force  vector  at  the  n'th  node  (the 
o 

point  of  measurement)  in  terms  of  the  displacement  vector 
at  the  O'th  node  (the  point  of  impact)  from  equation  (B.l) 

-  o;  «„  -  w0 

where 

•  •••  nK»,i2+  o*]''[K,],  +  i ,  O*]'1  (C.5 

In  our  case,  all  the  sections  are  identical  and  all  of  the 
K,,,  matrices  are  equal. 

C.2  An  alternative  method  of  computing  matrix  [R^]w 

Instead  of  having  to  invert  the  matrices  at  each  node  as 
shown  in  equation  (C.5)  we  could  successively  substitute 
from  the  left  hand  side  of  equation  (C.l)  into  the  right 
side  at  each  node.  This  would  give  the  forces  and 
displacements  at  the  O'th  node  in  terms  of  those  at  the 


n'th  node: 


which  is  only  one  inversion.  Unfortunately,  at  large  and 
,  the  hyperbolic  terms  in  matrices  G,  C2  and  F  , 
appendix  E,  are  very  large  and  almost  equal.  As  a  result 
the  elements  of  (Tw,  +  TW2  Q*  )  are  the  small  differences 
of  large  numbers  and  the  matrix  is  very  ill-conditioned. 
Even  though  this  is  only  a  2x2  matrix  for  the  symmetric 
case  we  lost  more  than  10  digits  for  three  multiplications 
at  frequencies  as  low  as  3000  Hz.  As  a  result,  the 
successive  inversion  method  was  used  for  all  cantilevers. 
Even  so,  some  care  must  be  taken  since  we  lose  6  digits  at 


10000  Hz. 


APPENDIX  D 


Derivation  of  the  Matrix  S 

The  equations  for  the  incident  waves  traveling  towards  the  right 
are 


Vj  =  exp  [  j  (wt-  £x)J 
u_j=  U , exp  j  (wt-  V  x)] 


(D.l) 


where  u  and  v  are  the  deflection  in  the  positive  direction  of  x 
and  y  respectively  (figure  2.5.2). 


For  the  reflected  wave  traveling  towards  the  left 
vR  =  v2  exp  [j(ut+  pxj]  +  V3  exp  [px+  j  to  t) 
uR  =  U2  exp  [  j  (wt+  r  x)] 


(D.2) 


At  node  zero  the  transmitted  displacements  are  v#  ,  and  u0 
which  are  related  to  the  internal  forces  at  the  node  by 


N  = 


W# 


Equation  (2.6.4)  written  for  its  all  components  is 


-jv,  =  jv2  +  V3  +  e0  /p 


(D.3) 


U,  -  u2  +  u0 


and  that  for  the  equations  of  equilibrium, 


Mx/Eip2  =  MR/EIp2  +  M,/EIp* 
Ft/EIp3  =  FR/Eip3  +  Fo/EIp3 
Pj/yEA  =  PR  /  y  EA  +Pfl/fEA 


(D.4) 


The  relations  between  the  components  of  NR  and  those  of 


displacement  W.  are 

Kr 

'  M,/Eipl  -  -V,  I 
Nj  =  ~Fj /Eip3  =  -jV,  \ 
Pj/j'  EA  =  -jU, 


*  nr  = 


Mg/EIf  - 


F^/EIp^  = 
Pr/YEA  = 


-v2+  v3 

-  jV.  +  V5 


=  jo. 


(D.5) 


No  = 


M0/EIfn  fvt  1 

-F#/Eip3  V  -  Q*  )  e„/pl  -  Q*  w 


P,/r  EA 


We  now  solve  for  V2  and  V3  in  terms  of  and  v0  and  URin  terms 
of  U,  and  U#from  equations  (D.3).  We  then  use  these  to  find 
the  elements  of  N  in  terms  of  V,  U,  and  U.0  in  equations 


' 


(D.5).  We  substitute  for  V(  and  U,  in  terms  of  Mj,  Fj  and  Pj 
from  equations  (D.5)  and  finally  a  set  of  linear  equations 
for  M1#  Fj.  and  Pr  in  terms  ofv,  0e  and  IL, 

NJ  -  s  we 

where 

H  H  (D . 6 ) 

i«CC3.o 


APPENDIX  E 


Differentiating  equation  (2.2.4)  we  have 


P  =  EA  u /  =  EA  y  [-EjCos  y  x  +  EjSinyxJ 

0  =  v7  =  A,^sinh|gx  +  A^cosh  p  x  -  A3|S  sin  x  +  k^pcos  p  x 
M  =  El  v*  =  EI^J  £  A  t  cosh  ^fix  +  A^sinh^sx  -  Aacos  p  x  -  A^sinpxj 
F  =-EI  v,/>=-Elp3^A(  sinh  p  x  +  A^cosh^x  +  A^sin^gx  -  A^cos^gxJ 


Substituting  x=0,  solving  for  E( ,  ,  Af,A2,  A3  and  A^  in 
terms  of  u,  v  and  their  derivatives  and  then  solving  for 
the  elements  of  Z;  at  x  =  L,  we  get: 


Zu  =  B  Z.Cc-0 
Zzi  =  G  ZzC^-0 


where 


B 


cos  o(| 


sin  o(, 


sin  <x, 


cos 


(E.l) 


*1  =  Cojh<Kz  4  CosVj 
<*1  =  +  Sm«fa 

Q 5  i  CojKW4  -  Ce&«<2 

=  Suihtfj  -SwMrfj. 


d,  ^2  S 


6*S 


az  flj 

A,  a2 

aj  aH  A| 


(E.2) 
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and 


o/.  =  —  ,  of3 


u  f  AW: 

=  '  p  =  Sr 


If  we  now  make  the  change  to  local  coordinates  for  the 
cantilevers,  and  put  in  the  boundary  conditions  for  the 
cantilever 


P=F=H=0atx2=h 


(figure  2.3.1) 


we  get  P2-/yEA  =  tano<5  u,* 


(E.3) 


Mu  /EI(SJ 


*Fa£  /EI  f 


=  A 


®|L  /p 


(E.  4 ) 


A=  J- 


COSka'14  Sin«<4  -Sink^q  Cos  <*4 


A  I  _  (tcsh^4Swi«<4  -  (sw-iktf.4  s^netfi*) 


,  coh 

-  - 

3  Cr 


■  ,  ^4  -  |3h  A  =  (  l-t  cosh«^4  c-os<*4  ) 


a)  Symmetric  section  (m  =  1) 


Substituting  from  equations  (E.l),  (E.2),  (E.3)  and  (E.4) 
into  equation  (A. 3)  we  finally  get  that 


*/*.*•»' .  ,  •  .*•  . 


-  ,  *,  *  _  •.  *  «  *  .*  *  " 


I 


0 


!  G  +  F  j 


where 


D.  =  -Vt 


2u)  J _ 

f‘r  AC 


-COSof, 


-sin  of. 


F  = 


o 

a.  a4 
b-  a, 


O 

o 

aa, 

to, 


0 

0 

au, 


a  =  . 


Coshrfq  Sinrfq  -  Sinh^  CoSaty 
1  ■+  Cosbwq  COS  i| 


0 

O 

a  a3 
b 


^  =  C  Cot  b  o<4  ju,  4  Sinbo(4  Cos*q  ) 


t>  - 


b)  Asymmetric  section  (  m  =  2) 

Substituting  equations  (E.l),  (E.2),  (E. 
equation  (A. 4)  we  get 


I  B  +  D,  •  C„  I 


P  Cr  4<m  <*3 

CO 


and  (E.4) 


where 


0 


0 


to  '  . 
d2  *  |9r/  2  A 


Ct  = 


0C*  2A  d4 


0  0 


(E  9) 


E  =  I  d  coseC,  d  sin*, 


(E.10) 


d  =  {sinho^  sinota)/A 


c)  Antisymmetric  section  (m  =  3) 


After  substituting  equations  (E.l),  (E.2),  (E.3)  and  (E.4) 
into  equation  (A. 5),  the  transfer  matrix  for  part  II  is 
obtained  as 


= 


B  +  D2  !  -  C 2 


!  G  +  i  F 


(E.ll) 


*.  \  \  *.  \ 


*  v-  v 
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